NEURAL CIRCUITS 



ORIGINAL RESEARCH ARTICLE 

published: 10 April 2013 
doi: 10. 3389/fncir.2013. 00064 




Theta-frequency resonance at the cerebellum input stage 
improves spike timing on the millisecond time-scale 

Daniela Gandolfi 1 - 2 , Paola Lombardo 1 , Jonathan Mapelli 2 - 3 , Sergio Solinas 3 and Egidio D'Angelo 1 - 3 * 

' Neurophysiology Unit, Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy 

2 Department of Biomedical, Metabolic and Neural Sciences, University of Modena and Reggio Emilia, Modena, Italy 

3 Brain Connectivity Center, IRCCS Fondazione Istituto Neurologico Nazionale C Mondino, Pavia, Italy 



Edited by: 

Chris I. De Zeeuw, Erasmus MC and 
Netherlands Institute for 
Neuroscience, Netherlands 

Reviewed by: 

Leonard Maler, University of Ottawa, 
Canada 

Attila Szucs, Hungarian Academy of 
Sciences Hungary 

"Correspondence: 

Egidio D'Angelo, Neurophysiology 
Unit, Department of Brain and 
Behavioral Sciences, University of 
Pavia, Via Forlanini 6, 27100 Pavia, 
Italy 

e-mail: dangelo@unipv.it 



The neuronal circuits of the brain are thought to use resonance and oscillations to improve 
communication over specific frequency bands (Llinas, 1988; Buzsaki, 2006). However, the 
properties and mechanism of these phenomena in brain circuits remain largely unknown. 
Here we show that, at the cerebellum input stage, the granular layer (GRL) generates 
its maximum response at 5-7 Hz both in vivo following tactile sensory stimulation of the 
whisker pad and in acute slices following mossy fiber bundle stimulation. The spatial analy- 
sis of GRL activity performed using voltage-sensitive dye (VSD) imaging revealed 5-7 Hz 
resonance covering large GRL areas. In single granule cells, resonance appeared as a reor- 
ganization of output spike bursts on the millisecond time-scale, such that the first spike 
occurred earlier and with higher temporal precision and the probability of spike genera- 
tion increased. Resonance was independent from circuit inhibition, as it persisted with 
little variation in the presence of the GABAa receptor blocker, gabazine. However, circuit 
inhibition reduced the resonance area more markedly at 7 Hz. Simulations with detailed 
computational models suggested that resonance depended on intrinsic granule cells ionic 
mechanisms: specifically, K S \ 0VJ (M-like) and KA currents acted as resonators and the per- 
sistent Na current and NMDA current acted as amplifiers. This form of resonance may 
play an important role for enhancing coherent spike emission from the GRL when theta- 
frequency bursts are transmitted by the cerebral cortex and peripheral sensory structures 
during sensory-motor processing, cognition, and learning. 
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INTRODUCTION 

Brain activity is characterized by complex temporal patterns, 
which often take the form of coherent oscillations (Buzsaki, 
2006). These are organized over defined frequency bands giving 
raise to the electroencephalographic rhythms. Important among 
these is the theta-band, which characterizes various functional 
states encompassing deep sleep, attentiveness, learning, and vol- 
untary movement. In the latter case, the commands generated 
by the motor cortex are associated with increased power in the 
theta-band (around 6-9 Hz in humans), which is then trans- 
mitted to cortical and subcortical centers (Gross et al., 2005; 
Schnitzler et al., 2006, 2009). Interestingly, these same frequen- 
cies are exploited by the thalamo-cortical circuits to elaborate 
motor commands and estimate kinematic parameters in behav- 
iors like whisking (Ahissar et al., 2000; Szwed et al., 2003; Kle- 
infeld et al., 2006; Zuo et al., 2011). But how do these sig- 
nals interfere with the activity of the receiving networks? One 
main hypothesis is that the receiving network is tuned over the 
transmission frequency band through resonance, being therefore 
more efficiently engaged when the transmitted frequency pat- 
tern is recognized (Llinas, 1988). This resonance could emerge 
both from membrane mechanisms based on specific ionic chan- 
nels and from the synaptic connectivity of the neuronal circuit 
involved. 



The cerebellum shows theta-band activity correlated with that 
of the premotor and motor areas during bimanual voluntary 
tasks in humans (Gross et al., 2005; Schnitzler et al., 2006, 2009). 
Moreover, single-neuron responses in the cerebellum are corre- 
lated with theta-frequency activity in the sensory-motor cortex 
of awake rats (Ros et al., 2009). Extracerebellar theta-frequency 
activity (Harish and Golomb, 2010) can be conveyed to the cere- 
bellum through two main pathways. The inferior olive, which also 
shows self-sustained theta-band oscillations (Marshall and Lang, 
2004), can generate theta-frequency spike bursts in climbing fibers 
(Mathy et al., 2009), and complex spikes in Purkinje cells (Welsh 
et al., 1995; Van Der Giessen et al., 2008). The granular layer 
(GRL) could also be entrained into theta-rhythms by extracere- 
bellar activity (Pellerin and Lamarre, 1997; Hartmann and Bower, 
1998; Ros etal, 2009). 

It has been reported that the interneuron inhibitory network 
made by Golgi cells shows theta-band resonant properties based on 
interneuronal connectivity through gap-junctions (Dugue et al., 
2009). Moreover, it was shown that neuronal elements of the GRL 
show intrinsic resonance in the theta-band, as both the granule 
and the Golgi cells respond maximally to input currents at around 
6 Hz (D'Angelo et al., 2001; Solinas et al., 2007a,b). However, it 
remained unclear whether the GRL resonates to extracerebellar 
inputs and which mechanisms would be involved. Here we show 
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that patterned sensory stimulation of the whisker pad causes GRL 
resonance at 5-7 Hz in the rat cerebellum GRL with marginal 
involvement of the inhibitory network. Resonance enhanced spike 
generation in granule cells raising time-precision on the millisec- 
ond time-scale. Theta-band resonance in the GRL could play an 
important role to tune the cerebellar circuit over one of the main 
frequency band of cerebro-cortical activity during generation of 
voluntary movement and other central functional states. 

MATERIALS AND METHODS 

Recordings were performed in vitro and in vivo and mathematical 
simulations were done using computational models. Experiments 
were performed using Wistar rats at postnatal day P20-P24 [inter- 
nal breeding, Harlan (Indianapolis, IN, USA)] both in acute slices 
and in vivo under anesthesia. All experiments were conducted 
in accordance with international guidelines from the European 
Community Council Directive 86/609/EEC on the ethical use of 
animals. 

STIMULATION PATTERNS 

In order to investigate the frequency - dependent properties of 
the GRL response, stimulus trials were repeated at frequencies 
between 1 and 10 Hz (1Hz steps). Each trial was composed of 
30 bursts and was repeated in a pseudo-random order (10 times 
in vivo, 20 times in both VSD recordings, whole-cell record- 
ings and simulations) to prevent potential effects due to short 
or long-term adaptation processes and to reduce the effect of 
response variability. The stimuli in vitro consisted of mossy 
fiber bursts composed of 3 pulses at 300 Hz intra-burst fre- 
quency. The stimuli in vivo consisted of 30 ms air puffs, which 
are known to generate short high-frequency bursts of similar fre- 
quency in the mossy fibers (Chadderton et al., 2004; Rancz et al., 
2007). 

RECORDINGS IN ACUTE CEREBELLAR SLICES 

Patch-clamp and VSD imaging recordings were obtained from 
parasagittal cerebellar slices (D'Angelo et al, 1995; D'Errico et al., 
2009). Briefly, the rats were deeply anesthetized with halothane 
( Sigma- Aldrich, Saint Louis, MO, USA) and decapitated. The 
cerebellum was removed, the vermis isolated and fixed on the 
vibroslicer stage (Dosaka, Kyoto, Japan) with cyano-acrylic glue. 
Acute 220 u.m thick slices were cut in cold cutting solution contain- 
ing (in mM): 130K-gluconate, 15 KC1, 0.2 EGTA, 20 HEPES, and 
10 glucose, pH adjusted at 7.4 with NaOH. Slices were incubated 
at 32°C in oxygenated extracellular Krebs solution containing (in 
mM): 120 NaCl, 2 KC1, 1.2 MgS0 4 , 26 NaHC0 3 , 1.2 KH 2 P0 4 , 2 
CaCh, 11 glucose (pH 7.4 when equilibrated with 95% O2 and 
5% CO2), at least lh before recordings. In some experiments 
the solution was implemented with the GABAa receptor blocker 
20 (iM gabazine (SR9553 1 ; Sigma- Aldrich) , which was maintained 
throughout the recording session. Slices were transferred to the 
recording chamber on the stage of an upright microscope (Zeiss 
Axioskop F2S, Oberkochen, Germany) and perfused at 1.5 ml 
min" 1 with oxygenated Krebs solution at 32°C with a feed-back 
Peltier device (TC-324B, Warner Instruments Corp., Hamden, CT, 
USA). Slices were immobilized with a nylon mesh attached to a 
platinum ST2-wire. 



Patch-clamp recordings 

Whole-cell current-clamp recordings were made in whole-cell 
patch-clamp configuration from granule cells as reported previ- 
ously (D'Angelo et al., 1995; D'Errico et al., 2009). Recordings were 
obtained using Multiclamp 700B amplifier (Molecular Devices, 
Union City, CA, USA) (3-dB; cut-off frequency = 10 kHz). 
Recordings were subsequently digitized at 20 kHz using pClamp 
9 (Molecular Devices) and a Digidata 1322A A/D converter 
(Molecular Devices). Patch pipettes were pulled from borosili- 
cate glass capillaries (Hilgenberg, Malsfeld, Germany) and filled 
with the following solution (in mM): 126K-gluconate, 4 NaCl, 
15 glucose, 5 HEPES, 1 MgS0 4 * 7 H 2 0, 0.1 BAPTA-free, 0.05 
BAPTA-Ca 2 +, 3 ATP, 100 |im GTP; pH was adjusted to 7.2 with 
KOH. This solution maintained resting free-[Ca 2+ ] at 100 nM 
and had a resistance of 7-10 MQ, before seal formation. The 
stability of patch-clamp recordings can be influenced by mod- 
ifications of series resistance and neurotransmitter release. To 
ensure that series resistance remained stable during the record- 
ings, passive cellular parameters were extracted in voltage-clamp 
by analyzing current relaxation induced by a 10 mV step from 
a holding potential of — 70 mV. According to previous reports 
(D'Angelo et al, 1993, 1995; Silver et al, 1996), the transients 
were reliably fitted with a monoexponential function yielding 
membrane capacitance (C m ) of 3.9 ± 0.2 pF (n = 18), membrane 
resistance (R m ) of 2.2 ± 0.3 Gfi (n = 18) and series resistance (_R S ) 
of 18.0 ± 0.9 MQ (n = 18). The 3-dB cell plus electrode cut-off fre- 
quency was fyc = (27tJ? s cm) 1 = 2.3 ± 0.1 kHz (n=18) and did 
not significantly change during recordings attesting their stabil- 
ity. Mossy fibers were stimulated with a bipolar tungsten electrode 
(Clark Instruments, Pangbourne, UK) via a stimulus isolation unit 
and stimulation intensity (±4—8 V; 100 u,s) was raised until the 
EPSPs generated spikes in 10-50% of the responses at 1 Hz from a 
membrane potential between —70 and —60 mV (mean 64.2 ± 2.7 
n= 18). From a comparison with previous data (D'Angelo et al., 
1995; Sola et al., 2004), between one and three mossy fibers were 
stimulated per granule cell. 

VSD imaging 

The procedures for VSD imaging were the same as reported in pre- 
vious papers (Mapelli and D'Angelo, 2007; Mapelli et al., 2010a,b). 
The slices were incubated for 30 min in oxygenated Krebs solu- 
tion added with 3% Di-4-ANEPPS stock solution mixed with 
50% fetal Bovine Serum (Molecular Probes). The dye (Di-4- 
ANEPPS, Molecular Probes) was dissolved and stocked in Krebs 
with 50% ethanol (SIGMA) and 5% Cremophor EL (a Castor oil 
derivative, SIGMA). Perfusion of standard extracellular solution 
(2-3 ml/min) maintained at 32°C with a feed-back temperature 
controller (Thermostat HC2, Multi Channel Systems, Reutlingen, 
Germany) was performed during the recording session. In a series 
of experiments, the extracellular solution was implemented with 
the GABAa receptor blocker, 20 |xM gabazine (SR95531; Sigma- 
Aldrich), which was maintained throughout the recording session. 
Mossy fibers were stimulated with a bipolar tungsten electrode 
(Clark Instruments, Pangbourne, UK) via a stimulus isolation unit 
using an electronic stimulator (STG 1008, Multi channel systems). 
The recording chamber was installed on an upright epifiuores- 
cence microscope (BX51WI, Olympus, Europa Gmbh, Hamburg, 
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Germany), equipped with a 20X (XLUM Plan FL 0.95 NA) (see 
Tominaga et al., 2000). The light generated by a halogen lamp 
(150W, MHF-G150LR, MORITEX Corp., Tokyo, Japan) was con- 
trolled by an electronic shutter (modelO, Copal, Co., Tokyo, Japan) 
and then passed through an excitation filter (X = 530± 10 nm), 
projected onto a dichroic mirror (X = 565nm), and reflected 
toward the objective lens to illuminate the specimen. Fluorescence 
generated by the tissue was transmitted through an absorption 
filter (X> 590nm) to the CCD camera (MICAM Ultima, Sci- 
media, Brainvision, Tokyo, Japan). The whole imaging system 
was connected through an I/O interface (Brainvision) to a PC 
controlling illumination, stimulation, and data acquisition. The 
final pixel size was 5 [im with 20 x objective. Full-frame image 
acquisition was performed at 1 kHz. The signal-to-noise ratio was 
improved by averaging 30 consecutive sweeps at the stimulus repe- 
tition frequency of 0.2 Hz. Given maximal AF/Fq ~ 1% and noise 
SEM ~ ± 0.1% (n = 12 slices), the signal-to-noise (S/N) ratio was 
about 10 times ensuring a reliable measurement of peak response 
amplitude. Data were acquired and displayed by Brainvision soft- 
ware and signals were analyzed using routines written in MATLAB 
(Mathworks, Natick, USA). The analysis was performed by evalu- 
ating the maximum response of the VSD signal obtained at each 
tested frequency (1-10 Hz) normalized to the first peak response. 

RECORDINGS IN VIVO 

Extracellular field recordings were performed in the GRL of 
Crus-IIa in 20-24-days-old Wistar rats (internal breeding, Har- 
lan). The rats were deeply anesthetized with urethane (1.4 mg/kg; 
Sigma-Aldrich; see Roggeri et al, 2008) dissolved in 0.9% NaCl 
and injected intraperitoneally. The heart rate (360-420/min) and 
respiratory rate (100-120/min) of each animal were constantly 
monitored and remained stable throughout the experiments. Ani- 
mals were placed on a stereotaxic table (David Kopf Instruments, 
Tujunga, CA, USA) covered with a heating pad and the body 
temperature was maintained at 38 ± 0.5°C through a feed-back 
temperature controller (Fine Science Tools Inc.; Foster City, CA, 
USA). The exposure of the cerebellar surface was performed fol- 
lowing previously reported surgical procedures (Bower and Wool- 
ston, 1983; Morissette and Bower, 1996; Lu et al, 2005). Briefly, 
a craniotomy was made on the right hemisphere (AP —13 mm 
ML 3 mm). The dura mater was carefully removed and extracellu- 
lar Krebs solution was placed onto the cerebellar surface (Roggeri 
et al, 2008). 

Local field potential (LFP) recordings from the GRL were 
obtained at the depth of 500 |xm with glass borosilicate pipettes 
(Hilgenberg) filled with NaCl 2M (0.5-1 Mfi). Insertion of the 
electrodes was at a 45° angle. Sensory stimulation was performed 
through a plastic pipette connected to a MPPI-2 Pressure Injec- 
tor (Applied Scientific Instrumentation, Eugene, OR, USA) and 
positioned 2-3 mm from the whisker pad. The low-frequency 
stimulation protocol delivered 30 ms air puffs (40psi) at 0.1 Hz 
frequency. Extracellular currents were recorded with Multiclamp 
700 A amplifier (Molecular Devices). The signals were band-pass 
filtered between 100 Hz (high-pass) and 2kHz (low-pass), digi- 
tized at 50 kHz using a Digidata 1322A interface, and stored on 
a PC using Clampex 9 (Molecular Devices). The same board and 
software were used to monitor and record body temperature and 



heart and respiratory rate and to generate stimulation pulses. The 
GABAa receptor blocker, 20 [iM gabazine (SR95531,Tocris Cook- 
son), was dissolved in Krebs solution, placed in a microliter syringe 
(Hamilton, Bonadus, Switzerland) and injected 500 [im deep into 
Crus-IIa of the cerebellar cortex close to recording electrode. 
Extracellular signals were acquired and processed off-line with the 
Clampfit 9.2 software (Molecular Devices). The signal-to-noise 
ratio was improved by averaging 100 consecutive traces and by dig- 
ital filtering (Gaussian low-pass filter 1 kHz). The LFPs consisted 
of two main waves, T and C. T was generated by rapid afferences 
through the trigeminal-cerebellar pathway, C by delayed affer- 
ences passing through the thalamo-cortico-ponto-cerebellar path- 
way (Morissette and Bower, 1996). In order to achieve a precise 
assessment of the LFP amplitude, only T measures were consid- 
ered for resonance analysis. LFP peak amplitudes were measured 
relatively to the preceding 10 ms baseline. 

COMPUTATIONAL MODELING 

Computational modeling was performed using a reduced version 
of a large-scale cerebellar GRL network model (Solinas et al., 20 10) . 
The granule cell (D'Angelo et al, 2001) and Golgi cell (Solinas 
et al., 2007a) models contain an explicit representation of ionic 
conductance mechanisms, which have been extensively tested with 
respect to large sets of experimental data. In order to generate the 
several seconds of neuronal activity required for these simulation, 
the modeling strategy was to use a single granule cell and a sin- 
gle Golgi cell in different combinations. The number of Golgi 
cells impinging on the granule cell (0-4) was simply simulated by 
linearly scaling the synaptic weight (see inset to Figure 6). The 
simulations were performed in the Python environment using 
models written in NEURON and were run on a 144 processors 
cluster (SiComputer INTEL MFSYS25V2). 

In brief, as in the previous network model (Solinas et al, 2010), 
the synaptic mechanisms included neurotransmission dynam- 
ics based on a vesicular release process generating short-term 
facilitation and depression. Release probability was derived from 
experimental works showing average values of 0.4 at the mossy 
fiber - granule cell synapse (Sola et al, 2004; D'Errico et al, 2009) 
and 0.4 at the Golgi cell - granule cell synapse (Mapelli et al, 
2009). Both the excitatory and inhibitory synapse were endowed 
with spillover mechanisms allowing activation of slow synaptic 
responses (Sargent et al., 2005). For the sake of simplicity, the 
synapses of the same kind (either excitatory or inhibitory) were 
set with identical release probability and postsynaptic receptor 
density and properties. In order to mimic synaptic noise and gen- 
erate jitter in spike timing, release at the mossy fiber - granule cell 
synapse was set in its stochastic version (Arleo et al, 2010) and 
background activity was generated in the mossy fibers activating 
the granule cell and the Golgi cell. 

The simulations were run by activating a single granule cell in 
all its 20 different synaptic combinations formed by four mossy 
fibers and four Golgi cell axon fibers (including the case of no active 
Golgi cell synapses to simulate inhibition block) (cf. Figure 8A). 
The Golgi cell was connected to the same four mossy fibers imping- 
ing of the granule cell plus other four, for a total of eight mossy 
fibers. All the mossy fibers carried low- frequency (4 Hz) back- 
ground activity. Then, in each simulation, 1-4 mossy fibers were 



Frontiers in Neural Circuits 



www.frontiersin.org 



April 2013 | Volume 7 | Article 64 | 3 



Gandolfi et a 



Theta-frequency resonance in the cerebellum 



activated by a 3 pulses - 300 Hz burst repeated at frequencies 
between 0.2 and 10 Hz (in 0.2 Hz steps) for 50 times (during the 
activation bursts, in these mossy fibers the background activity was 
interrupted). The data relative to membrane potential, membrane 
conductance, and current were stored and analyzed off-line using 
Matlab routines (Mathworks Inc.). 

The synaptic combinations determining the excitatory/inhibitory 
(E/I) balance of granule cells (Table 1) were considered to give 
specific relative contributions to the ensemble GRL response. The 
weight of these contributions was derived from estimates obtained 
in in vivo LFP recordings (Roggeri et al., 2008; Solinas et al, 2010; 
Diwakar et al., 2011) and from computations of the probability 
of connection between mossy fibers, Golgi cells, and granule cells 
(Solinas et al, 2010). These values are reported in Table 1. 

RESONANCE ANALYSIS 

In each experimental series (either LFP, VSD, patch-clamp, or 
modeling), for each tested frequency, the first five responses in 
each trial were excluded in order to allow the responses to sta- 
bilize. When required, the remaining 25 responses were averaged 
over all the trial repetitions at the same frequency. Resonance plots 
from LFP and VSD data were directly obtained by the response 
amplitude at the different frequencies. The response amplitudes 
at each frequency/; were normalized with respect to the value mea- 
sured at 1 Hz, y(fi) = {y(fi) — y(l Hz))/y(l Hz). Resonance plots 
for single-neuron responses (either whole-cell recordings or sim- 
ulations) were obtained by analyzing the frequency-dependent 
changes in spike count, first spike delay, first spike standard devi- 
ation, and average maximum depolarization (sc, sd, ssd, amd). In 
this case we calculated the Resonance Index (RI). At all tested fre- 
quencies, sc, sd, ssd, amd were measured, normalized between the 
extreme values for each cell and rescaled between 0 and 1. Thus, RI 
in a cell approaches one when the number of emitted spikes and 
membrane average depolarization are maximum and the delay 
and time-dispersion of the first spike are minimum. Then the four 
RI parameters (sc, sd, ssd, amd), were summed yielding a com- 
pound RI ranging from 0 to 4 and representative of resonance in 
the given cell. 

The resonance plots were usually well fitted by a double Lorentz 
equation y(f) (Siebert, 1986) (OriginPro8.0, Microcal Inc.) of the 
form: 





y(f) = 
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Table 1 | Cross-table of excitatory and inhibitory inputs to a g 
cell. 
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Exc/lnh (%) 


GoCO 


' GoC 1 


GoC 2 


GoC 3 


GoC 4 


Total 


MF 4 
MF 3 
MF 2 
MF 1 
Total 


0.43 
7.9 
0 
0 

8.33 


2.16 

12.78 

2.42 

0.97 

18.33 


3.97 

1.17 

24.84 

5.32 

35.3 


3.62 

1.48 

15 

6.29 

26.39 


I. 66 
0.51 
7.22 
2.26 

II. 65 


11.84 
23.84 
49.48 
14.84 
100 



The values indicate the relative contribution of each synaptic combination {from 

Solinas et al., 2010; Diwakar et al., 2011). 



Where, yo is a baseline level and, for each ith curve, A; is the 
underlying area, f c i is the peak frequency, and to,- is the width at 
y(f c i)/2. Lorentzian fittings allowed to find the resonance frequen- 
cies (-R-Fi_2 =/ c i,2) and to determine the relative weight of the 
two component of resonance plots with respect to their ampli- 
tudes y(RF2)/y(RFi). Lorentzian fittings allowed also to calculate 
RA = y(f c ) tw, which corresponds to the quality factor Q in reso- 
nance literature and characterizes a resonator's bandwidth relative 
to its center frequency. The goodness of fit was assessed with x 2 
statistics and by calculating the squared correlation factor, R 1 . In 
some cases, the identification of peaks did not require fittings so 
that RF and RA were determined directly from the raw data by 
using a peak-detection routine (when the two procedures were 
compared, the data obtained by peak-detection did not signifi- 
cantly differ from those obtained using Lorentzian fitting: e.g., see 
Figure 2B). 

STATISTICAL ANALYSIS 

Data are reported as mean ± SEM and compared for their statisti- 
cal significance by paired Student's f-test, unless stated otherwise 
(a difference was considered significant at p < 0.05). 

RESULTS 

RESONANT RESPONSES EVOKED IN THE CEREBELLAR CORTEX IN VIVO 

The response of the GRL to sensory inputs was evaluated by using 
repetitive air-puff stimulation of the whisker pad in urethane anes- 
thetized rats. A single air-puff is known to cause a short spike burst 
in mossy fibers (Chadderton et al, 2004; Rancz et al., 2007) and 
to determine the LFPs recorded from the GRL (Roggeri et al, 
2008; Diwakar et al., 2011) (Figure 1A). In order to investigate 
the frequency-dependence of LFPs in response to tactile stimu- 
lation, stimulus sequences of 30 impulses (30 ms duration) were 
delivered between 1 and 10 Hz in random order. After allowing 
responses to stabilize during the first 5 impulses, the amplitude 
of the average LFP obtained from the remaining 25 impulses was 
measured (Figure 1A). In control recordings, the LFP increased at 
5-7 Hz and then decreased again toward 10 Hz. Thus, the response 
was resonant in the middle of the theta-band, with a peak at 
7 Hz (65.6 ± 26.2% increase compared to 1 Hz; p < 0.01, n = 9) 
(Figure IB, left). 

The resonant properties of the GRL could reside into the 
intrinsic properties of granule cells (D'Angelo et al., 2001) but 
could also be shaped by the inhibitory Golgi interneuron net- 
work (Solinas et al, 2007b; Dugue et al, 2009). In order to test 
the role of the inhibitory circuit, the GABAa receptor blocker, 
20 jxM gabazine, was microperfused into the GRL. After 5 min 
from commencing perfusion, the LFP increased by 49.4 ±8.1% 
(n = 9; data obtained at 0.1 Hz, not shown), as expected from 
blockage of inhibitory transmission onto granule cells (Roggeri 
et al., 2008). Then, the stimulus sequence at different frequen- 
cies was repeated. The response was resonant in the middle of the 
theta-band, with a peak at 7 Hz (123.9 ± 4.9% increase compared 
to 1 Hz; p < 0.01, n = 9) (Figure IB, right). Apparently, resonance 
increased with respect to control and became more pronounced 
at 7 Hz compared to 5 Hz. 

Fittings to resonance plots were performed using double 
Lorentzian functions (see Eq.l) under the hypothesis that 
two components were present in the data distributions and 
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control 

1 Hz 7 Hz 10 Hz 




Freq(Hz) 

FIGURE 1 | Granular layer (GRL) resonance induced by sensory 
stimulation. (A). Schematic representation of the recording set-up in vivo. 
The lower inset illustrates that brief air-puff stimuli repeated at different 
frequencies are expected to generate repeated bursts in the mossy fibers 

(Chadderton et al., 2004; Rancz and Hausser, 2006; Rancz et al., 2007). By 
changing the frequency of sensory stimuli the frequency of the bursts is also 
changed. The GRL response is revealed by local field potential (LFP) 
recordings. The plot shows the time-course of LFP amplitude in response to 



gabazine 



1 Hz 7 Hz 10 Hz 




1 23456789 10 

Freq(Hz) 

subsequent pulses at three different stimulus frequencies. (B) LFPs recorded 
from the GRL at different frequencies of sensory stimulation in control and 
after intracerebellar microperfusion of gabazine. The normalized maximum 
LFP amplitude is reported in the plots. Note resonance at 5-7 Hz, which is 
maintained in the presence of gabazine. Data are reported as mean±SEM. 
Red lines fitting the data point represent double Lorentzian functions (broken 
lines are the two individual component functions) obtained on control and 
gabazine data, respectively. 



that their relative contribution could explain shape changes 
caused by gabazine application. In control, fittings yielded 
1^1 =4.4 ±0.9 and RF 2 = 6.7 ±0.3 with y(RF 2 /y(RFi) = 2.2 
(R 2 = 0.86; x 2 = 0.02). In the presence of gabazine, fittings yielded 
£f 1 =4.7±1.6 and RF 2 = 6.7 ±03 with y(RF 2 )/y(RFi) = 4.8 
(_R 2 = 0.9; x 2 = 0.003). Thus, gabazine application modulated the 
resonance profile by raising the relative weight of the component 
peaking around 7 Hz with respect to that peaking around 5 Hz. 

RESONANT RESPONSES EVOKED IN ACUTE CEREBELLAR SLICES 

Voltage-sensitive dye (VSD) imaging was performed in parasagit- 
tal cerebellar slices in order to define the spatial distribution of res- 
onance in the GRL. Mossy fiber stimulation was performed using 
3 pulse-300 Hz bursts repeated at frequencies varying between 1 
and 10 Hz in random order. The signals generated by mossy fiber 



stimulation invaded the GRL (Figure 2A) reproducing patterns 
reported previously (Mapelli et al., 2010a,b). In the same slices, 
20 |iM gabazine application increased the response to mossy fiber 
stimuli (average AFAFrj increase of 42.7 ± 11%; p < 0.01 n = 4), as 
expected from the disinhibitory effect of the drug (Mapelli et al, 
2010a,b). 

The GRL responses showed maximum AF/Fnat input fre- 
quencies of 5-7 Hz (Figure 2A). This behavior was quanti- 
fied by constructing resonance curves in ROIs covering limited 
areas of the GRL (or even single pixels), which showed char- 
acteristic profiles peaking in the same 5-7 Hz frequency range 
(Figure 2B). In the presence of gabazine, resonance at 5-7 Hz was 
still evident. 

Resonance maps were constructed by plotting the resonance 
frequency (RF) measured in each pixel on color-scale (Figure 3A). 
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FIGURE 2 | The spatial distribution of resonance in the granular layer 
(GRL). (A) Optical maps of VSD responses in the GRL evoked by repeated 
bursts stimulation delivered to the mossy fibers bundle (mf) at three different 
frequencies. The VSD maps were generated by taking the maximum 
depolarization for each pixel both in control and after gabazine perfusion. 



Scale bar 100|.im. (B)The panels show the resonance curves corresponding 
to specific ROIs of the GRL indicated in (A) (squares). The data points show 
maximum responses either at 5 Hz (green) or 7 Hz (orange), both in control 
and after gabazine perfusion (arrowheads indicate points at the resonance 
frequency). Broken lines are single Lorentzian functions fitted to the points. 



The maps clearly showed that the large majority of single pixel 
resonance frequencies were in the 5-7 Hz region. After gabazine 
application, RF appeared still concentrated at 5-7 Hz, although 
the area covered by 7 Hz RF became prevalent. 

The average pattern of changes was evaluated by gen- 
erating average resonance curves, which were constructed 
from all single pixel resonance curves in each of 11 slices 
(Figure 3B). The average resonance curve showed two sub- 
peaks at 5 Hz (38.6 ±9.3% increase vs. 1 Hz) and 7 Hz 
(31.1 ±7.2% increase vs. 1 Hz), which were both signifi- 
cantly higher than nearest neighbor points (p < 0.05, paired 
(■-test). In the presence of gabazine, the average resonance 
curve increased compared to controls (137.1 ±9.6% increase; 
p > 0.4 f-test paired with nearest neighbors), with a maxi- 
mum at 7 Hz (Figure 3B, right). In control, fittings yielded 
£Fi = 5.0±0.6 and RF 2 = 7.0±0.6 with y(RF 2 )ly(RFi) = 0.7 
(R 2 = 0.92; x 2 = 0.04). In the presence of gabazine, fittings yielded 
£Fi=4.4±0.5 and RF 2 = 7.0 ± 0.4 with y(RF 2 )ly(RF x ) = lA 
(R 2 = 0.99;x 2 = 0.006). The ensemble changes caused by gabazine 
in the resonance curve were to increase the absolute size of the 
component peaking around 7 Hz and its relative weight with 
respect to that peaking around 5 Hz : similar to what observed 
in the LFP in vivo. 

The extension of 5-7 Hz resonance in the GRL was assessed by 
counting the number of pixels with the same resonant frequency 



(Figure 3C). On average, 5 Hz resonance occurred in 23.4 ± 4.7% 
of the active area and 7 Hz resonance occurred in 13.4 ±2.9% 
of the active area. Overall, the only significant peak was that at 
5 Hz (p < 0.01; paired f-test with nearest neighbors), suggesting 
that this was the most represented resonant frequency. In the 
presence of gabazine, the amount of area showing resonance at 
7 Hz was increased compared to that at 5 Hz with a single signif- 
icant peak at 7 Hz (31.7 ± 5.8% of active pixels; p < 0.01; paired 
f-test with nearest neighbors). In order to determine the origin 
of these differences, we measured resonance in individual pix- 
els (evaluated using resonance amplitude, RA: see Materials and 
Methods), which showed similar amplitude in control and after 
gabazine perfusion (p < 0.00005, comparison of 50 pixels in con- 
trol and after gabazine chosen at random). Therefore, since single 
pixel resonance was not significantly affected by gabazine (see 
above), the amplification of the resonance curve with gabazine 
was probably due to the increased area over which resonance 
occurred. 

In aggregate, these results show that the cerebellar GRL in acute 
slices generates maximum responses when the input bursts con- 
veyed through the mossy fibers are organized in theta-frequency 
range. The results closely resemble those obtained from LFPs 
in vivo, in that theta-frequency resonance persists but tends 
to increase at 7 Hz compared to 5 Hz after blocking inhibitory 
synaptic transmission. 
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FIGURE 3 | Synaptic inhibition modulates granular layer (GRL) 
resonance. (A) The maps show the resonance frequency (RF) for each of 
GRL pixels before and after application of gabazine. The maps reveal a 
prevalence of areas showing resonance at 5 and 7 Hz and an increase of the 
7 Hz area after gabazine. Scale bar 100 |im. (B)The resonance profiles 
obtained by all active GRL pixels and averaged over several recordings (n = 11 
control and n = 4 gabazine) reveal statistically significant peaks in the 5-7 Hz 
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range. The amplitude of the resonance plot increases after gabazine. 
(C) Histograms report the relative areas covered by each RF averaged over 
several recordings (n= 11 control and n = 4 gabazine). Note that the main 
peak in control occurs at 5 Hz but moves to 7 Hz in the presence of gabazine. 
Data are reported as mean ± SEM. Red lines fitting the data point represent 
double Lorentzian functions (broken lines are the two individual component 
functions) obtained on control and gabazine data, respectively. 



RESONANT RESPONSES EVOKED IN SINGLE GRANULE CELLS BY 
PATTERNED SYNAPTIC ACTIVITY 

In order to investigate the cellular basis of resonance in the GRL, we 
performed whole-cell recordings from granule cells in acute cere- 
bellar slices. The granule cells, which showed high input resistance 
(_R in = of 2.2±0.3Gfi «=18) and low membrane capacitance 
(C m = 3.9 ± 0.2 pF, n=18), generated repeated spike discharge 
upon depolarizing current injection (D'Angelo et al., 1995, 2001) 
(Figure 4A, inset). The granule cells were maintained at a holding 
potential between —60 and —70 mV, from which mossy fiber low- 
frequency stimulation could elicit spikes in <50% of responses. 
Then, high-frequency burst patterns were repeated between 1 
and 10 Hz in random order and the modifications in synaptic 
excitation were analyzed (Figure 4A; left). 

In most recordings (n = 7), at 5-7 Hz, the burst depolariza- 
tion increased along with the probability of generating spikes, 



which also became more precise and occurred with shorter delay 
(Figure 4A; left). In order to obtain an estimate of resonance 
from these measurements we calculated the RI (see Materials and 
Methods for definition) through the combination of the following 
parameters: spikes count (sc), first spike delay (sd), first spike SD 
(ssd), and average maximum depolarization (amd). It should be 
noted that these parameters are strongly correlated, so a stronger 
depolarization commonly elicited stronger spiking and shorter 
latency of the First spike (absolute values for sc and ssd are reported 
in Figure 4A, inset). Interestingly, Rl/frequency plots for all the 
parameters (sc, sd, ssd, amd) showed a resonant shape and com- 
monly peaked between 5 and 7 Hz (Figure 4B; left). The RIs were 
then summed to obtain a compound RI for each cell (Figure 4C, 
left). A similar analysis was performed in recordings in which 
gabazine was applied (n = 7). The granule cells became more 
excitable and made more spikes generating protracted discharges, 
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FIGURE 4 | Resonance in spike emission from single granule cells. 

(A) Left; Single sweeps from a single granule cell recorded in control 
solution are shown over imposed at three different frequencies (2, 5, 
9 Hz). Note the increased number, shorter delay and higher precision of 
spikes at 5 Hz. Blue traces show the average depolarization. Note that 
maximum depolarization is obtained at 5 Hz. Right; Single sweeps from 
a single granule cell recorded in the presence of gabazine are shown 
over imposed at three different frequencies (2, 5, 9 Hz). Note the 
increased number, shorter delay, and higher precision of spikes at 5 Hz. 



Blue traces show the average depolarization. Note that maximum 
depolarization is obtained at 5 Hz. The inset in the middle shows granule 
cell responses elicited by 1 50 ms-2pA current steps. The inset on the 
right shows the spike count and the spike variance in control and after 
gabazine perfusion. (B) Left; In a single granule cell, the Resonance 
Index (Rl) for individual parameters in control (sc, ssd, sd, amd) revea 
peaks between 5 and 7 Hz. Right; In the presence of gabazine Rl reveal 
peaks between 5 and 7 Hz (C) Compound Rl in the same granule cell in 
control (left) and in the presence of gabazine (right). 



but the resonance patterns remained unvaried. In most cases, the 
granule cells generated more spikes and the first spike occurred 
earlier and with higher precision around 5-7 Hz (Figure 4A; right) 
causing a higher maximum average depolarization at the same fre- 
quencies (Figure 4A; right). Accordingly, the Rl/frequency plots 
showed main peaks at 5-7 Hz both for individual parameters (sc, 
sd, ssd, amd) (Figure 4B; right) and for the compound Rl index 
(Figure 4C; right). 



The probability of observing a given RF for each of the para- 
meters (sc, sd, ssd, amd) is reported in Figure 5A, which shows 
that the maximum concentration of resonance frequencies occurs 
at 5-7 Hz, both in control and in the presence of gabazine. When 
the probability of occurrence were summed, a clear resonance 
curve peaking at 5-7 Hz was generated. Another way to repre- 
sent the behavior of the granule cell population was to average 
their compound RIs. Also this average compound Rl/frequency 
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FIGURE 5 | The global effect of resonance on spike emission in the 
granular layer network. (A) The histogram shows, for each frequency, 
the probability of observing a certain RF for spike number, spike delay, 
spike standard deviation and maximum average depolarization (n = 7). 
The superimposed plots show, at each frequency, the sum of 
observations in individual parameters. Note peaks at 5-7 Hz. (B)The plot 



shows the compound resonance index (Rl) calculated from all the 
neurons recorded in control and in the presence of gabazine. Note peaks 
at 5-7 Hz. Data are reported as mean ± SEM. Red lines fitting the data 
point represent double Lorentzian functions (broken lines are the two 
individual component functions) obtained on control and gabazine data, 
respectively. 



curve showed major peaks at 5-7 Hz (Figure 5B). In control, fit- 
tings to these plots yielded: RFi = 4.8 ± 0.3 and RF 2 = 7.0 ± 0.2 
vnthy(KF 2 Vy(RF\) = 0.8 (.R 2 = 0.93; j 2 = 0.06). In the presence 
of gabazine, fittings yielded RF\ = 2.7 ± 0.5 and RF 2 = 7.4 ± 0.6 
with y(RF 2 )/y(RF 1 ) = 1.1 (£ 2 = 0.86; x 2 = 0.1). Thus, the reso- 
nance profile of the ensemble responses obtained from several sin- 
gle granule cells showed that, after the application of gabazine, the 
relative weight of the component peaking around 7 Hz increased 
with respect to that peaking around 5 Hz. 

In summary, the ensemble of spike-related parameters in gran- 
ule cells generated resonance compatible with that observed in 
global network measurements obtained with LFPs (cf. Figure 1) 
and VSD imaging (cf. Figures 2-3), suggesting that resonance can 
be traced to elementary cellular phenomena mostly residing in the 
granule cells. 

MODELING GRANULAR LAYER RESONANCE 

The mechanisms underlying circuit resonance were explored by 
performing simulations with a model derived from a previous 
large-scale network of the GRL (Solinas et al., 2010). The model 
included detailed realistic representations of granule cells and 



Golgi cells and the related synapses (D'Angelo et al, 2001; Nieus 
et al, 2006; Solinas et al, 2007a,b; Arleo et al., 2010). Similar to 
what observed experimentally (D'Angelo et al., 1995; Sola et al, 
2004), the mossy fiber - granule cell synapses were endowed with 
stochastic neurotransmission mechanisms, so that simulated gran- 
ule cell responses were affected by noise and showed jitter in spike 
generation (Figure 6A). Once subjected to input patterns identi- 
cal to those used in whole-cell recordings, the granule cells in the 
model generated Rl/frequency plots for the resonance parameters 
(sc, sd, ssd, amd) and for compound RI, with main peaks in the 
5-7 Hz region (Figure 6B), both in control and when simulat- 
ing the "gabazine" condition (phasic and tonic inhibition blocked; 
Figure 6C). 

The granule cells in the model received all the basic combi- 
nations of excitatory and inhibitory connections occurring in the 
real network, amounting to just 20 different combinations (see 
Table 1). Indeed, an active granule cells receives 1-4 excitatory 
and 0-4 inhibitory synapses, on average (Solinas et al., 2010). 
The Rl/frequency plots of resonance parameters (sc, sd, ssd, amd) 
for all the synaptic combinations are shown in Figure 7A. These 
plots reveal that, due to the different excitatory/inhibitory (E/I) 
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FIGURE 6 | Modeling resonance in spike emission from single granule 
cells. Resonance in a model granule cell embedded into the granular layer 
network model. The traces and plots on the left show the case of two 
excitatory and two inhibitory synapses and simulate a control experimental 
condition. The traces and plots on the right have two excitatory synapses, 
while tonic and phasic inhibition are blocked simulating experimental gabazine 
application. The mossy fibers in the model are activated with patterns 
homologous to those used for investigating resonance in acute cerebellar 
slices (3-spike bursts at 300 Hz repeated 50 times at frequencies ranging 
from 1 to 10 Hz in 0.5 Hz increments). The inset shows a schematic of the 



model, with afferent fibers carrying background only (gray) or background 
interrupted by the bursts (black). (A) The last 45 of 50 sweeps from a single 
granule cell in the network model are shown superimposed at three different 
frequencies. Note the increased number, shorter delay and higher precision of 
spikes at 6 Hz (control and inhibition blocked) compared to the other two 
frequencies. The average traces are shown superimposed in color. (B) In the 
same model granule cell as in (A), the R.I. for individual parameters (left: sc, 
ssd, sd; right: amd) reveal peaks in the theta range. (C) Compound R.I. 
obtained from the data shown in (B). Both in (A,B) the peak around 7 Hz is 
enhanced after blocking synaptic inhibition. 



balance, individual granule cells do not show identical resonant 
properties, although the Rl/frequency plots usually peak in the 
5-7 Hz range. Interestingly, resonance was more marked (evalu- 
ated using resonance amplitude, RA: see Materials and Methods) 
and more precisely centered at 5-7 Hz for E/I < 1 (Figure 7B). 
Moreover, when synaptic inhibition was blocked (either in the 
transient or both in the transient and tonic component) , RF shifted 
toward 7 Hz. These observations provide a possible explanation 
for the variability in resonance properties of single granule cells 
measured in whole-cell recordings and suggest that GRL reso- 
nance emerges from the statistical distribution of microscopic 
parameters characterizing granule cells with different E/I. 



In order to evaluate how the different E/I balance in granule cells 
influenced RF and Rl/frequency plots, we determined the statisti- 
cal distribution of these parameters (Figure 8A). The probability 
of observing a given RF for each of the resonance parameters (sc, 
sd, ssd, amd) is reported in Figure 8A for two different probabil- 
ity distributions. A uniform distribution, in which the probability 
of finding the different E/I balances is identical, was compared 
to the probability distribution reported by Diwakar et al. (2011) 
from the analysis of LFP recordings in vivo and VSD recordings in 
acute cerebellar slices. In both cases, the highest concentration of 
resonance parameters was observed at 5-7 Hz, with a higher peak 
using the Diwakar's distribution. 
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FIGURE 7 | Contribution of different E/l combinations to 
single-neuron resonance. The impact of the excitatory/inhibitory balance 
on resonance is shown for all E/l combinations, indicated as £„//„ (n is the 
number of synapses impinging on a given granule cell). Simulations were 
carried out as in Figure 7. (A) The effect of E/l combinations on Rl for 
individual parameters (sc, ssd, sd, amd) are shown in the four panels. The 
columns / 0 correspond to absence of phasic inhibition. Note that peaks in 
the plots more often occur between 5 and 7 Hz. (B)The plots reports RA 



(defined in Materials and Methods) for individual parameters (sc, ssd, sd, 
amd) in the cases of coexistence of E and / (left), of phasic inhibition 
blocked (middle) or of both phasic and tonic inhibition blocked (right). For 
each point, RF is reported using a color code. Note that 5-7 Hz RF clusters 
are at E/l < 1 and correspond to the highest RA values, consistent with the 
enhancing effect of inhibition on resonance. Note also that, when both 
phasic and tonic inhibition are blocked, RF shifts toward higher values, as 
observed experimentally. 



The behavior of the whole granule cell population was 
reconstructed by summating and averaging the compound RIs 
for a given statistical distribution of E/I combinations. The 
Rl/frequency curve showed major peaks at 5-7 Hz, and the 
7 Hz peak became dominant after blocking synaptic inhibition 
(Figure 8B). The reason of the enhanced effect of inhibition at 
7 Hz compared to 5 Hz was analyzed by applying RI analysis to 
Golgi cells in the simulated network. In response to mossy fiber 
bursts, Golgi cells had a main resonance peakat around 6-7 Hz, 
thereby depressing the granule cell response more effectively at 
this specific frequency (Figure 8B). 

MODEL PREDICTIONS OF THE CELLULAR MECHANISMS OF 
RESONANCE 

Intrinsic granule cell resonance can be elicited by sinusoidal cur- 
rent injection depending on an M-like current (K s i ow ) and and 
A-current (KA) acting as "resonators" and a persistent Na cur- 
rent acting as "amplifier" (Hutcheon and Yarom, 2000; D'Angelo 
et al., 2001). In the present simulations, during repetitive burst 
transmission, -K" s i ow was higher at high-frequency while KA 



was higher at low- frequency, so that the two current-frequency 
curves intersected and generated a minimum ,K" s i ow + KA cur- 
rent at 5-7 Hz. At the same frequency, the persistent Na cur- 
rents (Nap) and the NMDA current showed enhanced acti- 
vation. The combined effect was that of generating a surplus 
of 2pA in the 5-7 Hz range (Figure 9A). It should be noted 
that, over the high input resistance of granule cells (in the 
Giga-ohm range), a 2pA current is indeed capable of caus- 
ing a substantial enhancement in membrane depolarization and 
spike generation, as shown in previous studies (D'Angelo et al, 
2001). 

The time-course of ^ s i ow and KA currents is shown in 
Figure 9B. While -K" s i ow accumulated as the frequency of stimu- 
lation increased toward 10 Hz, KA decreased toward 10 Hz. The 
summed current was maximal at 5-7 Hz, just at the time of 
granule cell spike burst generation. The mechanism of frequency- 
dependence of -K"si ow and KA is analyzed in Figure 9B. K s \ ow 
activated slowly with time constants in the 100 ms range and 
therefore benefited of short inter-burst intervals. Conversely, KA 
inactivated rapidly with time constants in the 10 ms range and then 
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FIGURE 8 | Modeling the global effect of resonance on spike emission 
in the granular layer (GRL) network. The global effect of resonance on 
spike emission in the GRL network has been modeled using different 
proportional contributions of granule cells with respect to their Ell 
balance. (A) The histograms show, for each frequency, the probability of 
observing a certain RF for individual parameters (sc, ssd, sd, amd). The 
superimposed plots show, at each frequency, the sum of observations in 
individual parameters. Note peaks at 5-7 Hz as in experimental recordings. 
In the histogram on the left we assumed a uniform probability distribution 



of Ell combinations. In the plot on the right we assumed the probability 
Ell distribution derived from Diwakar et al. (2011 ) in vivo. (B)The plot 
shows the compound Rl calculated from model simulations using the Ell 
distribution derived from Diwakar et al. (2011) either including all E/l 
combinations (black trace) or with phasic and tonic inhibition blocked (gray 
trace; gabazine effect). Note the similarity of Rl plots to those observed 
experimentally. The compound Rl of a Golgi cell in the network (dashed 
trace) shows a major peak at 6-7 Hz, explaining why inhibition block 
enhances Rl at the same frequencies. 



took time to de-inactivate, so that KA benefited of long inter-burst 
intervals. The gating properties of -K" s i ow and KA were therefore 
mechanistically correlated with generation of resonance in granule 
cells during mossy fiber burst transmission. 

DISCUSSION 

Following recognition of the role of oscillations and resonance 
as fundamental aspects of neuronal communication in the brain 
(Llinas, 1988; Buzsaki, 2006), their phenomenological proper- 
ties and mechanisms have remained unclear in several circuits. 
This paper demonstrates that the cerebellum GRL, once acti- 
vated with periodic inputs, shows resonance at 5-7 Hz. Reso- 
nance was manifest both in vivo following sensory stimulation 
with air puffs delivered to the whisker pad and in acute cere- 
bellar slices following electrical stimulation with short bursts 
delivered to the mossy fiber bundle. Resonance was modified 
but persisted after blocking inhibitory synaptic transmission. 
Thus, theta-frequency resonance largely depended on intrin- 
sic properties of the neurons involved (Hutcheon and Yarom, 
2000). Interestingly, resonance was expressed by a change in 
the composition of spike bursts emitted by granule cells with 



a modulation of spike timing on the millisecond time-scale. 
Therefore, GRL resonance could recode with millisecond pre- 
cision the theta-bursts at the input into new bursts at the 
output. 

GENERAL PROPERTIES OF GRANULAR LAYER RESONANCE: THE 
RELATIONSHIP WITH CIRCUIT INHIBITION 

A striking property of GRL resonance evoked by theta-frequency 
patterns is that it has almost identical RF and gabazine sensitivity 
in vivo and in vitro. Resonance in vivo and in vitro are mecha- 
nistically correlated since each sensory stimulus generates a short 
high-frequency burst in mossy fibers (Chadderton et al, 2004; 
Rancz et al., 2007), so that the theta-frequency sensory volley 
is mimicked by theta-frequency stimulation of the mossy fiber 
bursts (Roggeri et al., 2008). This suggests that resonance observed 
in vivo is almost entirely generated in the GRL circuit rather than 
in the up-stream sensory pathway passing through the trigeminal 
nucleus (De Zeeuw et al, 1996). Also cerebro-cortical compo- 
nents were excluded, since the present analysis concerned only the 
trigeminal T wave of the cerebellar LFP (Morissette and Bower, 
1996). 
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FIGURE 9 | Ionic mechanisms controlling single-neuron resonance. The 

ionic mechanisms controlling single granule cell resonance in response to 
patterned synaptic activity were analyzed in simulations on the E 2 /'i 
combination. (A) The values of currents generated by /C st)w , KA, Nap, and 
NMDA channels in response to simulated granule cells bursts were measured 
20 ms after burst activation at three selected frequencies (2, 5, and 9 Hz). 
While / K A decreases, / /e st)w increases with frequency, so that i w + i K slow has 
a minimum at around 5 Hz. The ;' Na and / NMnA increase at 5 Hz, so that a net 
inward current of about 2 pA is generated at 5 Hz. (B) The plot on top shows 



/ K a and i K slow during simulated granule cells bursts delivered at frequencies 
of 2, 5, and 9 Hz. Note the enhancement of the total current at 5 Hz. The plots 
in the middle show the kinetics of K alow activation and KA inactivation. Note 
that K stow accumulated activation with frequency, whereas KA accumulates 
inactivation with frequency, explain the global resonant effect of their 
combination. The bottom tracings show membrane potential at the three 
different frequencies. All traces are taken from the simulations with two 
active excitatory and inhibitory synapses (average of 25 of 30 consecutive 
traces, the first 5 excluded). 



Both in LFPs in vivo and in VSD and whole-cell record- 
ings in slices, the application of gabazine, a GABA-A receptor 
antagonist, to block inhibitory transmission between Golgi cells 
and granule cells did not suppress resonance but rather enhanced 
its component at 7 Hz. In general, GRL resonance was appro- 
priately described by two components peaking at 5 and 7 Hz, 
and the 7 Hz peak prevailed when inhibitory transmission was 
blocked. In simulations, the Golgi cell showed maximum activity 
at 6-7 Hz, suggesting that the inhibitory circuit was indeed more 
effective at this particular frequency. It is thus possible that areas 
with 7 Hz resonance correspond to circuit regions that are less 
strongly inhibited than those with 5 Hz resonance. Double reso- 
nance peaks have been reported in the inferior olive, where the 
relative weight of the higher frequency peak is also depressed by 
synaptic inhibition (Llinas and Yarom, 1986). It could therefore 
be that in these circuits the inhibitory circuit has a modula- 
tory role on resonance occurring along the main transmission 
pathway. Double resonance was also observed in cortical neu- 
rons, but in that case one peak was in the theta and one in the 
gamma band (Cobb et al, 1995). The inhibitory circuit passing 



through the cerebellar glomerulus also involves GABA-B receptors. 
Although a potential role of GABA-B receptors in GRL resonance 
remains unknown, it should be noted that GABA-B receptors 
down-regulate GABA-A receptor-mediated transmission at the 
Golgi cell - granule cell synapse [both by reducing GABA release 
(Mapelli et al, 2009) and postsynaptic receptor activation (Bran- 
dalise et al., 2012)] and enhance granule cell intrinsic excitability 
[by reducing an inward rectifier current (Rossi et al., 2006)]. At 
present, a major contribution of GABA-B receptors to resonance 
is therefore improbable. 

GRANULAR LAYER RESONANCE EMERGES FROM MILLISECOND 
CONTROL OF SPIKE TIMING 

Single cell analysis revealed that resonance in LFP and VSD record- 
ings reflects the microscopic nature of spike generation in gran- 
ule cells. At the RF, the granule cells showed higher probability 
of emitting spikes, which occurred with higher precision and 
shorter delay, resulting in a larger average depolarization. Since 
the LFP in the GRL is mostly due to extracellular spike cur- 
rents (Diwakar et al., 2011), a higher probability of generating 
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spikes and a higher synchrony of the spikes could well explain LFP 
resonance. As far as the VSD signal is concerned, since this tech- 
nique is currently unable to precisely measure individual spikes, 
VSD resonance was more probably correlated with average cell 
depolarization. 

PREDICTION OF RESONANCE MECHANISMS USING REALISTIC 
COMPUTATIONAL MODELING 

The mechanisms of granule cell resonance during synaptic trans- 
mission were explored using a realistic computational models 
of the GRL circuit. The model, faithfully reproduced all macro- 
scopic properties of resonance, including RF at 5-7 Hz and the 
shift toward 7 Hz caused by inhibition blockage. At the micro- 
scopic level, the model showed that resonance properties depended 
from the combination of single cell components in a statistically 
distribute manner and that the specific combination of exci- 
tatory and inhibitory synapses activating granule cells in vivo 
was especially effective in generating resonance in ensemble net- 
work activities. In the model, intrinsic granule cell resonance 
emerged from the interplay of K s \ ow and KA causing a mini- 
mum repolarizing current at 5-7 Hz. However, opposite to the 
case of sinusoidal current injection (D'Angelo et al, 2001), -K" s i ow 
prevailed at higher frequencies while KA prevailed at lower fre- 
quencies. This happened because, different from sinusoidal cur- 
rents, the bursts have fixed duration and what determines the 
frequency-dependence of channel activation and inactivation is 
the inter-burst interval. Thus, in subsequent bursts, -K" s i ow (Hu 
et al, 2002) accumulated activation while KA accumulated inac- 
tivation. Granule cell resonance was amplified by the persistent 
Na current (Hu et al., 2002) and by voltage-dependent unblock 
of the NMDA current in the just-subthreshold region, which 
enhanced EPSP-spike coupling (D'Angelo et al., 1995; Mapelli and 
D'Angelo, 2007). As already pointed out (D'Angelo et al, 2001) It 
should be noted that, due to the high input resistance of granule 
cells, currents of just a few picoamperes in the threshold region 
resulted in significant effects on membrane depolarization and 
spike generation. 

RESONANCE AND OSCILLATIONS IN THE GRANULAR LAYER 

Resonance is a condition occurring when a physical system under- 
goes a periodic activation with a frequency equal or close to the 
intrinsic oscillation frequency of the system itself, so that such a 
system tends to oscillate at its maximum amplitude. The nature 
of oscillations and resonance depends on the physical details of 
the system involved. So, what is the relationship between reso- 
nance and oscillations in the GRL? Theta-frequency oscillations 
are observed in the GRL during resting activity in the awake rat 
and monkey (Pellerin and Lamarre, 1997; Hartmann and Bower, 
1998; Courtemanche et al, 2009). Computational analysis indi- 
cates that these oscillations require an intact feed-back inhibitory 
loop (Maex and De Schutter, 1998; Solinas et al, 2010). Conversely, 
as we report here, GRL resonance reflects intrinsic properties of 
granule cells. Finally, resonance and oscillations in the inhibitory 
interneuron network require gap-junctions between Golgi cells 
and other network components (Forti et al., 2006; Dugue et al., 
2009) and possibly also intrinsic pacemaking in Golgi cells and 



other network components (Forti et al, 2006; Galliano et al, 2013). 
Thus, the circuit appears to be composed of two sub-systems, a 
resonator (the mossy fiber - granule cell synapse) coupled to a 
resonant oscillator (the Golgi cell inhibitory network). This lat- 
ter also provides synchronicity through lateral inhibition and can 
enhance resonance in the granule cell population. In aggregate, 
resonance can amplify the granule cell output when the mossy 
fiber input is conveyed at 5-7 Hz. At this frequency the inhibitory 
circuit can spontaneously oscillate, thereby creating a condition 
at which the system can optimize phase locking and information 
transmission. 

As pointed out above, by demonstrating similar resonance 
properties in vivo and in vitro, our results suggest that the major 
resonance mechanisms revealed by whisker pad stimulation reside 
in the GRL network. Indeed, the mossy fiber bursts evoked by eye- 
movements (Kase et al, 1980) correlate strictly with the sensory 
stimulus without any dependence on the phase or duration of the 
stimulus pattern. Nonetheless, the neuronal networks providing 
inputs to the GRL and even the whiskers might also present forms 
of resonance. For example, the thalamo-cortical circuit operates 
on a low-frequency bandwidth during whisking (Ahissar et al., 
2000; Szwed et al., 2003; Kleinfeld et al, 2006) and both vibrissal 
motoneurons (Harish and Golomb, 2010) and trigeminal neurons 
(Wu et al., 2001) show forms of low- frequency oscillations and 
resonance. Therefore, low-frequency resonance in the cerebellar 
circuit is probably part of a resonant system distributed over the 
different nodes of the sensorimotor network controlling whisk- 
ing. It should also be noted that mechanical resonance of vibrissae 
occurs at high-frequency and is probably relevant for detect- 
ing sharp changes in objects surface like edges or irregularities 
(Hartmann et al, 2003). 

CONCLUSIONS AND FUNCTIONAL IMPLICATIONS 

Granular layer network resonance, by reflecting millisecond regu- 
lation in spike emission, could fully exploit the outstanding timing 
capabilities of granule cells (Cathala et al, 2003; D'Angelo and De 
Zeeuw, 2009; Diwakar et al., 2009, 2011) and fine-tune informa- 
tion transfer (Arleo et al., 2010). One potential implication of 
GRL resonance is that, by occurring on the same band of Purk- 
inje cell and inferior olive oscillations (Llinas and Yarom, 1986; 
Welsh et al, 1995; Lang et al, 2006; Abrams et al, 2012), could 
help maintaining a high level of coherence in the activity of the 
whole olivo-cerebellar system. On a larger scale, the cerebellum 
could be optimally designed to detect information carried by 
cerebro-cortical theta cycles implementing a well-tuned transmit- 
ter - receiver system (D'Angelo et al., 2009). For example, motor 
commands for whisking are emitted at theta-rhythm entraining 
the thalamo-cortical circuit (Szwed et al, 2003) and the cerebel- 
lum (O'Connor et al, 2002; Lang et al, 2006) providing the basis 
for complex resonant loops. Since oscillations in the GRL and in 
the Purinje cell - inferior olive circuits are tunable between 5 and 
20 Hz (i.e., beyond the theta-band), the 5-7 Hz resonance reported 
here maybe restricted to specific aspects of behavior. A hint comes 
from the specific sensitivity of the GRL to theta-bursts conveyed 
through sensory and cortical inputs, which can induce long-term 
synaptic plasticity in the cerebellum (Roggeri et al, 2008; Diwakar 
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et al, 201 1 ) as well as in the cortex and hippocampus (Larson and 
Lynch, 1988; Huerta and Lisman, 1995). An intriguing hypothesis 
is that resonance may improve learning of salient patterns in rela- 
tion to voluntary movement (Gross et al., 2005; Schnitzler et al, 
2006), attention and sleep, in which theta activity prevails. 
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